SWENDSEN-WANG DYNAMICS ON 7L d 
FOR DISORDERED NON FERROMAGNETIC SYSTEMS. 
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Abstract. We study the Swendsen-Wang dynamics for disordered non ferromagnetic Ising models on 
cubic subsets of the hypercubic lattice Z d and we show that for all small values of the temperature 
parameter f3~ 1 the dynamics has a slow relaxation to equilibrium (it is torpid). Looking into this 
dynamics from the point of view of the Markov chains theory we can prove that the spectral radius goes 
to one when the size of the system goes to infinity. This means that, if we want to use the Swendsen- 
Wang dynamics for a computer simulation, we have a slow convergence to the stationary measure in low 
temperature. Also it is a good example of a non-local dynamics that relaxes slowly to the equilibrium 
measure. 
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1. Introduction 

The study of the Swendsen-Wang (SW) dynamics has received a considerable attention in recent 



years; see the original paper [E6fl and also |15l [19, 20 , [7j , where some rigorous mathematical results on 



the convergence speed to the equilibrium measure have been obtained. 

The Swendsen-Wang algorithm was proposed for the simulation on a computer of the Ising or Potts 
model (not necessarily ferromagnetic) to obtain a more rapid simulation with respect to the standard 
Monte Carlo methods. In this paper we are precisely interested in the behavior of the algorithm for non 
ferromagnetic Ising models defined on cubic subsets of 7L d . We do not deal with general Potts model, but 
we conjecture that many of the results proved here are true, in some sense, also in that case. 

The SW dynamics is a reversible Markov chain having the Gibbs Ising measure as the stationary one. 

A relevant feature of the SW algorithm is that it is non-local, in the sense that many spins can change 
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value simultaneously. This feature allows us to hope for a fast convergence to equilibrium. In fact for 
ferromagnetic systems with an external field on Z d it is proved (see |2(| ) that this method of simulation 
has a fast approach to the equilibrium measure; it is effective also for small values of the temperature 
parameter where standard Monte Carlo methods are typically slow. On the other hand, it has been 
recently proved that the speed of convergence to the equilibrium measure of the SW algorithm for the 
mean field is not always rapid . Recently it was proved that the mixing time of the SW dynamics for 
the Potts model in a box of Z d of side length L, at the critical temperature, with a large enough number 
of states, goes to infinity as L — > oo; see where some explicit inequalities are established for the mixing 
time of the SW algorithm and for the mixing time of the more standard Glauber dynamics. 

We will prove that the SW dynamics for disordered non ferromagnetic Ising systems, on 1 d , spin-glass 
systems, has a slow approach to the equilibrium measure for large values of the parameter (3 (the inverse 
of the temperature). More precisely we prove that the spectral gap of the associated Markov chain 
decreases to zero when the side length of the box goes to infinity if (3 is large enough. 

Some arguments and results of the paper can be framed in the dynamical Griffiths singularity. The 
dynamical Griffiths singularity for the ordinary Glauber dynamics has been studied by various authors 
(see U H [| g [[§ Notice also that our results have some similarities with [0, in the sense that we 
show another model (but with a large region of the parameter temperature in which we can apply the 
proofs) where the SW algorithm has a slow convergence to equilibrium for the physically relevant case of 
the lattice Z d . 

We notice that in |l8| it is proved that the Swendsen-Wang dynamics for ferromagnetic Ising model in 
the low temperature case is non-ergodic and hence there is a slow convergence to the equilibrium in finite 
volume. In any case there are no results on the spectral gap and we conjecture that the spectral gap 
remains bounded away from zero as the side length goes to infinity for the SW dynamics applied to the 
ferromagnetic Ising model. Thus we believe that it is not easy to replace our construction to prove that 
the spectral gap shrinks to zero. Furthermore, we conjecture that the the behavior is not governed by 
the occurrence of large ferromagnetic blocks, almost isolated, as it is in some other cases (see for example 
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H ||, ||). Also for these reasons we believe that this kind of Griffiths singularity is a new phenomenon 
not studied yet. 

Our results are true in any dimension for all the Ising models where the interactions are random 
variables and belongs to the interior of the support. This hypothesis on the support of the measure 
can seem unnatural but there are a lot of probability measures that respect those conditions; for example 
the uniform measure on an interval that includes the origin or the Gaussian measures, so the results are 
highly applicable. 

We stress that the SW dynamics is related to the Fortuin-Kasteleyn (FK) representation that has an 
important role to study the phase transition. In fact for ferromagnetic models there is a 1-1 relationship 
between bond percolation and spin phase transition (see |lj, |l)) but for frustrated systems there is a 
weaker relationship that allows us only to establish inequalities on the phase transition for the spin-glass 



(see y, |S|, O, |14, 21 ). In any case it is not known for the general spin-glass if there is a phase transition or 
not. Thus our results are not proved using arguments on the equilibrium measures or the phase transition 
of these Ising models, but we use directly the structure of the SW dynamics. 



2. Swendsen-Wang dynamics and main definitions 

In this section we will recall the definition of the Swendsen-Wang dynamics and the main definitions 
that will be used in the paper. Our goal is to prove that the spectral radius associated with the Markov 
chain representing the SW dynamics goes to one when the size of the system grows to infinity invading 
all the lattice 1 d . In order to show the essential line of the proof, we will not find explicit bounds on 
the spectral radius but we will only prove that it goes to one. In any case all the estimations could be 
done explicitly. So, if one is interested in usable bounds for the spectral radius, it is possible to perform 
explicit calculations following our proofs. 



Let us consider the lattice Z d where the edges are all the couples of vertices x, y having ||a; — y\\ = 1 
and || • || is the Euclidean distance. 
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Some sets. Given a set of edges Bc£wc define the vertex set of B, V(B), as the set of all vertices 
v G V such that there exists w G V with {v,w} G B. Given a set of vertices A we define E(A) as 
the set of edges e e E such that V^({e}) G A; we remark that given a set of vertices V\ C V we have 
V(.E(Vi)) C Vi and given a set of edges Ei C E the relation £?(V(£i)) = E\ is true. In the following if 
we do a set operation using a vertex set A and an edge set B in the same time, then we will see these 
two sets as vertex sets; for example if A C V and B C E then An B corresponds to An V(B). 

Boundaries. For a set of vertices V C Z d let us denote by dV, the boundary ofV, as the set of vertices 
belonging to V such that for every vertex u G V there is a vertex w G V c such that {w, w} G E is an 
edge. 

For a set of vertices V we define also the edge boundary dbV as follows e = {it, u} G dbV if and only if 
e n V ^ and e n V" c # 0. 

o o o —A J) i— 

OOOO — ©OOG— 

oooo — eooo— 

O O O (j) |— 

(a) (b) 

Figure 1 . The boundary of a set of vertices in (a) and its edge boundary in (b) 



We will deal with interactions J = {Jb}beE(z d ) that are i-i-d. random variables with a distribution Qj 
having a neighborhood of zero inside its support (supp Qj). 

Our proof works for all (3 G (f3h, oo) where flh depends on the distribution Qj and on the dimension d 
of the lattice. We could also generalize the proof to some other kinds of distributions of the interactions, 
for example the uniform measure on (—6, —a) U (a, b) for suitable a, b > 0, but only for some limited 
intervals of the parameter [3 and not for all (3 G (f3h, oo). 

To fix the ideas and the notation we can think that Qj is a continuous distribution whose support 
includes [— 1, 1]. If this is not the case we can re-scale the support of Qj such that [—1, 1] C supp Qj; 
this will correspond to change the value of the parameter [3. 
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We will consider an (hyper) cube of vertices Bu{v) with length of the side Al and center v. Let us fix 
an event in the realization of interaction on the cube Bn{v) having positive probability. 

For the edges e G dbB 4 i we require that J2eed b B it \^ e \ < V^ 2 ; f° r every edge e G E(B 4 i) \ dbB 2 i we 
set J e G (1 — l/l d , 1]; finally for every e G <9;,i?2Z we define the interactions through these rules (only for I 
large enough): 

1. We choose d faces Si, S 2 , ■ ■ ■ , Sd among the 2d faces of the cube B21 in such a way that every face 
Si is not the image of another of these faces using R, the central reflection of the cube. The faces Si 
are sets of vertices. The central reflection of the cube R can be used as a function on the vertices 
of the cube into itself or as a function on the edges of the cube into itself (it is an isomorphism of 
the graph). So for every edge e = {vi, v 2 } of the cube also Re := {Rvi,Rv 2 } is an edge. 

2. For every edge e G dbB 2 i having a vertex v = (vi, v 2 , ■ ■ ■ , va) belonging to a face {Si}i = i....,d of the 
cube we give the following rule: if | X^=i v i\ 1& an even number we impose that J e G (ad — l/l d , ad) 
otherwise Jb G (—sad, —sad + l/^ d ) where s is in (0, 1) (in the following we will specify it) and ad 
is an opportune positive constant. 

3. For every edge e G dbB 2 i that we have not already defined we use this rule: if Jn e G (a^ — l/l d , ad) 
then J e G (-sad, -sa d + l/l d ) and if J Re G (-sad, -sa d + l/l d ) then J e G (a d - l/l d , a d ). 

So all the interactions inside the cube B41 are defined. 

We have fixed the interactions on a finite region so that there is a positive probability to find a cube 
B41 of this kind. When are fixed all the interactions on the graph Z d one can find a cube B41 of this 
kind at a finite distance from the origin Qj-almost surely (by ergodicity with respect to the translation). 
We will indicate with the box [-L, L] d n Z d . Now, using standard arguments, it is easy to see that 
there exists a subsequence of boxes with L(l) G N strictly increasing with / such that only for a 

finite number of I the box does not contain a cube B41 with the desired configuration of interactions 
(Qj-almost surely). 
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Configuration Space. Let us define the space of the spin configurations Q = { — 1, l} zd ; let us write 
Qa = { — 1, 1} A where A is a finite subset of Z d . We define the space H — {0, l}- 6 ^) (configurations of 
edges) and H\ = {0, 1} E< - A '> (configurations of edges on A). 

Ising model. Let us define the Hamiltonian of the Ising model on a finite subset A C Z d as: 

(1) H\ t j(a) = - ^2 Ji,kO-i<Jk- 

i,k:\\i-k\\ = l 

Here i, k are vertices in A n Z d and {i, k} G E(A), a — {cri}i e \ n z d is a configuration belonging to 
The Ising model (or Ising measure) on a finite subset A C Z d is defined as 

(2) k a , j0 (<t) = 

where Za,jp is a normalizing factor {partition function). Sometimes we will omit the dependence from 
the box A and the parameters J. 

Note that the interactions J = {J e } e eE(A) are random variables; in this case the system is called 
disordered or quenched since there will be a random measure or a stochastic process parameterized by 
some random variables. 

Swendsen- Wang dynamics. The Swendsen-Wang dynamics (or process) is a Markov chain that updates 
alternatively the configuration of edges Ha and the configuration of spins f^A. The SW dynamics has 
as space of states Q A x Ha with a finite A (we can consider a cube in d dimensions as A) and it is an 
algorithm generally used to simulate the Ising model on a finite space. 

The Markov transition probability Pp(u),r];u),rj), with uo.uo e Oa and 77,77 G Ha, that realizes the SW 
dynamics is not homogeneous, indeed in every odd time i e N the algorithm updates the configuration 
of the edges (in H A ) and in every even time ( £ H SW dynamics updates the configuration of the 
spins (in ^a)- We will also introduce the Markov transition probability M / | w/ (w, to), with w, cD e £1a, 
that is the Markov transition probability on the even times restricted to the space of the states J7a (so 
Mp W = Pp\n A )- For a configuration (cr, 77) e £!a x Ha at time t = In, with tieN, the algorithms update 
the configuration with these two steps: 
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SW1) This step updates the configuration of edges. Let A C E(A) be the subset of edges e = {x, y} £ E(A) 
such that J e a x a y > 0. For every edge e G A, independently of the others, with probability 
p e = 1 — e^l' 7 "! in the updated configuration rf e = 1 and with probability 1 — p e in the updated 
configuration 77^ = (we remark that if J e = then p e = 0) . 

SW2) This step updates the configuration of vertices. Let us consider the e-clusters, maximal connected 
component, of the graph where two vertices, u and w, are connected if there exists a path 7 = 
(u, ei, . . . , e n , u>) such that r\ ei — 1 for every i = 1, . . . , n. For every e-cluster C we will choose a 
vertex u and we will decide with probability 1/2 that o~ v is equal to 1 or — 1. For all the other 
vertices of the e-cluster the values of the spin are assigned following the rule J e a x a y > where 
e = {x, y}. This procedure define in every vertex the configuration Q\. 

It is immediate to see that the Markov chain M§ w is irreducible and aperiodic (since every element 
of the transition matrix Mp is larger than zero); np, the Gibbs measure of the Ising model, is the 
invariant measure of the Markov chain M^ w (see Edwards and Sokal |2(], |ll|). We will indicate with 

: Qa x Ha — > [0,1] the invariant measure of the Markov chain Pp. The measure vp is the FK-measure 
related to the Ising model (see f2l|j). 

It is easy to understand that also the transition probability Pp is irreducible and aperiodic on the 
states J7a x Ha. It is believed that the SW algorithm is faster than the standard Monte Carlo algorithm 
because this dynamics is not local so it is more easy to leave a metastablc state |2(| . 

Reversible property. The SW algorithm is a reversible Markov chain; it was proved in the original 
paper pfij ] . The reversibility or the detailed balance was directly used in ]2l| to show that the equilibrium 
measure of this algorithm is the Ising measure. So a consequence of reversibility and regularity (i.e. the 
Markov chain is irreducible and aperiodic) is that all the eigenvalues of the transition matrix are real and 
inside the interval (—1, 1]. These eigenvalues have multiplicity one with probability one if the interactions 
{ JejegE are absolutely continuous random variables. We order these eigenvalues 



(3) 



1 = A > Ai > • • • > A n > -1 
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where n is the number of configurations on a finite number of vertices V that is equal to 2' y L 

Exponential Time and the Spectral Radius. For the exposition of the exponential time we follow p3 . 
We will give the definitions for a system with a finite number of states. Let P be a regular transition 
matrix associated to a Markov chain and let ir be its unique invariant measure. Let / : J7a ->Kbea real 
function that is square-integrable with respect to the stationary measure ir , i.e. / € L 2 (ir). 

Let us write 

(4) n(f) = J2 /MttH 

wen 

for the mean value of / with respect to the measure 7r and 

(5) C ff (t)= ]T f(u)f{wMu)P t (»,&)-W)] 2 = 

(6) =VarM)~\ E (/M- = for, (/) -£(/,/) 

for the autocorrelation function. We are denoting with Varv(/) the variance of the random variable / 
and with £*(/, /) the Dirichlet form of the transition matrix P ; in the following the Dirichlet form will 
be used to bound the spectral radius of the transition matrix. We also notice that C//(0) = Var^{f). 
Let us define the normalized autocorrelation function as 

{<) PffW = r< (n y 

For finite regular Markov chains Pff(t) decays exponentially for large t; we define, for / 6 L 2 (7r),the 
exponential f- autocorrelation time 

(8) r expJ = lim sup — — — 

t^oo — log 

and the exponential autocorrelation time 

(9) 7~ eX p — SUp 'Texp.f- 

/GL 2 (7r) 

Thus, r exp is the typical time that the Markov chain employs to relax towards the equilibrium measure. 
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Equivalently one can see P as an operator on the Hilbert space L 2 (ir). Let define R the spectral 
radius of the operator P as 

(10) i? = sup|A 4 | 

If there is only the eigenvalue Ao = 1 on the unit circle then R for a finite space of states is strictly less 
than 1. The relation 

(11) R = e - 1/T ^ 

is in [ST89, Propositions 2.3-2.5] (see also (2^]). So if R is close to one the Markov chain relax slowly to 
the equilibrium measure. 

Our Markov chain is dependent from the size of the box Al = [— L,L] d n Z d , so the quantities 
T~exp,f — T exp,f,L, i~ exp = T exp ^, and R = Rl are function of L. We will choose some particular functions 



Jl to show, using the variational characterization, that Rl goes to one almost surely (see Theorem 4 
for L that goes to infinity (and hence T exp .L — ► oo). 

Magnetization and related events. Let us define the magnetization of a configuration a G fl\ on a finite 
region A C A of vertices as: 

(12) M(A,*) = ±-J2 a '» 

' ' veA 

so M(A, a) is a rational number belonging to the interval [—1,1]. 
Let us define the accordance on the set of edges B C E(A) as: 



(13) K (B,a) = J L 

Finally let's define the events: 



\B\ 

1 1 e={v,w}eB 



(14) S+(S) = {a : K(d b B 2l , a) >1-S}, 



(15) 



Sf(S) = {a : K{d b B 2l ,a) <-l + S}, 
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where S is a constant in (0, 1). 

It will not be difficult to see that for all S € (0, 1) and (3 large enough the probability of the event 
(5) tends to one when I goes to infinity. In Section ^ we will give some detailed estimations for the 

probability of the events 5 ; + (S) and S[~ (5) that will be used to prove that the spectral radius goes to one 

when A — » Z d . 

3. Satisfied edges and Peierls surfaces 
In this section we will give other definitions. 

Clusters. For every v € Z d let us denote by C v = C v (a) the cluster at the vertex v. C v is defined as 
the maximal connected set of Z d having constant sign. In this case two vertices u and w are connected if 
there is a path 7 — (u, ei, vi, . . . , u n -i ; e n , w) such that a u — o~ Vi and a u = a w for every i = 1, . . . ,n — 1. 
In an obvious way we can extend this definition of cluster to a generic graph Q = (V, E). 

Given a set of vertices V (resp. a set of edges E) we will call C a cluster on V (resp. on E) the cluster 
on the subgraph H = (V,E(V)) (resp. H = (V(E),E)). 

Peierls surface. We call a set of edges S a Peierls surface if there exists a configuration a such that 
S = dbC(a) : where C is a cluster. It is possible to associate a plaquette to every edge e of the lattice, 
where this plaquette is the (d — l)-dimcnsional unit cube orthogonal to e intersecting the edge in its 
center. It is a standard fact that the plaquettes associated to a Peierls surface form a closed (without 
boundary) (d — l)-manifold. 

Satisfied, unsatisfied and neutral edges. For a fixed configuration a we call satisfied an edge e = {u, v} 
such that J uv o~ u a v > 0. It is called unsatisfied if JuvO~ u o v < and neutral if J uv a u a v = 0. It is easy 
to see that given a set of vertices V if we consider a spin-flip on V (i.e. the sign of the spin of all the 
vertices in V changes) then all the edges e S dbV that were satisfied (resp. unsatisfied) become unsatisfied 
(resp. satisfied); all the other edges that are not in d^V or that are neutral remain of the same kind. Note 
that if J has an absolutely continuous distribution then the probability that there is a neutral edge is 
zero. 
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4. Spectral radius for the SW dynamics 



In the first two lemmas of this section we present some geometrical results. Lemma 4.1 is an easy 



calculation that shows an inclusion between two events. The following Lemmas 4.2, 4.3 are a consequence 



of some results on the minimal surfaces; we write the proofs following the reference In Lemmas 4.4 



4.6 and 4.7 we give the essential bounds on the probability that the Markov chain does a transition from 
any state in Sf(S) to any state in S^~(6) and a bound on the probability of the event Q\ (Sf(S) US^~(5)). 
In Theorem we combine the previous bounds and, using the variational characterization of the second 
eigenvalue of the transition matrix, we will show that the spectral radius goes to one so the exponential 
autocorrelation time goes to infinity. In all the bounds of the lemmas we do not use that the Markov 
chain is reversible but in Theorem [4.8 it is essential for the use of the variational characterization. 



Lemma 4.1. Let S £ (0, 1/2). If a £ Q \ (S^(5) U S l (5)) and I is large enough then at least one of the 
following inequalities is true: 

a) \M{dB 2U a)\ < 1 - J/4; 

b) \M(dB 2(l+lh a)\ < 1-6/ A. 



Proof. The proof of this Lemma is an algebraic calculation. Let us suppose that M(dB 2 i, a) > 0. If 
M (dB2i , er) < we can do a global spin flip of the configuration leaving unchanged the absolute values 
of the magnetization and the value of K. So we can consider only the case M(dB 2 i, cr) > 0. Now there 
are two different cases to analyze: 

a) M(dB mi) ,a)>0; 

b) M(dB 2{l+1) ,a) <0. 

First we prove the lemma in case a). As a consequence of our hypothesis K(dbB 2 i, u) < 1 — 6 so there 
are at least 8\dbB 2 i\/2 edges e = {v,w} on which the product of a v a w is equal to —1. So at least 
S\dB 2 i\/4 are the negative spins that are either on dB 2 i or on dB 2 ^ +1 y Therefore we have that either 
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M{dB 2l ,a) < 1-8/2 or 



M(dB 2(l+1) ,a) < ' 



dB 2 (i + i) 



\dB 2[l+1) \-2^\dB 2l \ 



It is clear that lirr^oo \dB 2 i\/\dB 2 ^ + i)\ = 1 so for ^ large enough \dB 2 i | / 1 c3i? 2 (z+i) I > 1/2 and we obtain 
the inequality of the lemma. Case b) follows by an analogous calculation just using the inequality 
K{d b B 2U a) > -1 + 5. n 

The following two lemmas are connected to problems of minimal surfaces for a given boundary. 

Lemma 4.2. Let S G (0, 1/2). If for a configuration a <E Q the inequality \M(dB 2 i, a)\ < 1 — S is verified 
then there exists a Peierls surface 7 inside box B 2 i where all the edges in 7 \ dbB 2 i are unsatisfied and 
|7 \ d^B 2 i \ > 5pd\dbB 2 i\, where pd is a constant in (0, 1] depending only on the dimension d. 

We will use also the following analogous lemma that we will not prove because the proof is basically 



the same of Lemma 4.2 



Lemma 4.3. Let 5 G (0,1/2). If for a configuration a 6 Q the inequality |M(9_B 2 (;+i), cr)| < 1 — 5 is 
verified then there exists a Peierls surface 7 inside B41 \ B 2 i where all the edges in 7 \ (dbB 2 i U dbBa) are 
unsatisfied and I7 \ (dbB 2 i U dbB±i)\ > 5pd\dbB 2 i \, where pd is a constant in (0, 1] depending only on the 
dimension d. 



Before the proof of Lemma 4.2 we give these two definitions 



(16) At (a) = {ve dB 2l : o v = 1}, 



(17) A-{o) = {»6 dB 2l : o v = -1} = dB 2l \ A+(a) 



We notice that, as an easy consequence of hypothesis of Lemma 4.2, there exists a constant p 
min{|, jL} g (0, -L] such that 



(18) 
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Proof of Lemma 4.2. To fix ideas we impose < M(dB2i,cr) < 1 — 5 but the other case — 1 + 5 < 
M(dB2i,cr) < is analogous. Let us consider the region B21 and the family of clusters C — {Q}i=i m 
on this region. The family of clusters C form a partition of the vertices in V(B2i); we fix an order on 
the clusters belonging to C. We will do a spin- flip of all the vertices belonging to the first cluster C G C 
which do not touch the boundary <9i?2z; sometimes in the following, doing this spin-flip of the vertices of 
a cluster C, we will say that we are eliminating the cluster C . We will fix a new order on this new family 
of clusters on B21 and we will continue to eliminate the first cluster that do not touch the boundary dB2i ■ 
In the end of this procedure we will obtain a configuration <7 (a) in which all the clusters on B21 touch the 
boundary dB 2 i- We will denote with = {Ci{(j { ^)} l=1 ... M(o ) the new family of clusters on i?2Z- It is 
easy to verify that the boundary of the clusters can only decrease with this procedure, so if we can find 



the Peierls surface 7 (verifying Lemma 4.2) in this new configuration this is allowed also in the old 
configuration a. So we will find a Peierls surface 7 such that: 

(19) 7CUAC,(a (o) )cUAC,M. 

Also it is easy to verify that Af(a) = Af(a {a) ), and so M(dB 2 i, a {a) ) = M(dB 2 i,a). Now we will 
analyze two sub-cases 

i) dC n dB, = 0, 

ii) dC n dB, ^ 0, 

where O is the center of the cube In the first case we consider 7 = dbCo and we do the projection from 
the point O of the surface 7 on <9i?2Z , we obtain the inequality of the Lemma with standard arguments 
(see |). 

Indeed in this first case we have: 

(20) ^<1- P 

\ObB2l\ 

in fact either dC C A+{<j^) or dC C A-(cr^). 

Before we analyze case ii) let us define the following partial order for the clusters on the cube B21: the 
cluster in the center of the cube Co is the smallest one (Co d C x for every x € B21) and two clusters 
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C x , C y are ordered C x < C y if for every path a starting in the origin and ending in a vertex of C y the 
path a has at least a vertex belonging to the cluster C x . 
If there is a cluster C such that 

1 j io s |as 2; | s 10 

we define 7 = dC, otherwise we start the following procedure. Let us choose a maximal cluster with 
respect to the partial order and consider the spin-flip of this cluster. We notice that with this procedure 
we can only decrease the boundary of the clusters inside the square B 2 i so after every such step up is 
more difficult to find a surface 7 verifying the Lemma; so if we can prove that there is such a kind of 
surface then there is the same surface also in the original configuration. Every time that we do a spin-flip 
of a cluster we will verify if there exist a cluster respecting the inequality (|2l]). When there is a cluster 
C verifying the inequality ( pl[ ) we will define 7 = dC. Surely there exists a time in which this procedure 
ends because the relation ( pi] ) can be not true only if for all the clusters 

f09 , \ddndB 2l \ P 

{ } \dB 2l \ < IO' 

So in the end of this procedure there is a unique cluster C that verify the inequalities 

P l^ n ^| ^ p 
[3) W < \dB 2l \ <5- 

Let us consider the intersection of the cluster C with dB m where m €{?,/ + 1, ... , 21}. Let us define 

(24) Pm := lC ^ B ! nl me {1,1 + 1, . . . , 21}; 

\OD m \ 

let's also write 

(25) ^:=inf— m, n €{/,/ + 1, 21}. 

m.n p n 

We remark that by definition v < 1. Again we will analyze two sub-cases 

a) 0<^< A ; 

b) &<!/<!. 
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The sub-case a) can be proved with the projection of the surface as in the case i) of the previous page 
(just changing the constant pj). 

We prove the sub-case b); on every level m G {/,..., 21} it is easy to see that jq < p m < 1 — jg- Now 
using the isoperimetric inequality on a compact manifold, see || pag. 209, for every level m € {I, . . . , 21} 
we obtain the inequality 



(26) 



\dCndB m \ > K d 



— \dB n 
10 1 



d-2 
d-1 



where K d is a positive constant depending on the manifold (in this case only on the dimension of the 



space because the manifolds are cubes). For dimension d — 2 the RHS of the inequality (26) becomes 
independent of dB m and it is immediate to understand that 



(27) 



\ddndB m \ > l. 



Therefore we have a lower bound for the measure of the surface S = (DC \ E(Bi)) C 7, S is the surface 
dC out of the square Bi (S is (d — 2) dimensional manifold but it is also a set of edges); so we can write 



(28) 



K d l inf 

me{l,...,2l} 



10 



\dB„ 



d-2 

" 1 7^ J—l 

> —K d pd^l 



d-l 



In the last inequality we use that p £ (0,1] so p*- 1 > p for every d > 2 and \dBi\ — 2dl d 1 . This 
concludes the proof for the last sub-case, so the lemma is proved for all d > 2. 

□ 

We fix the value of the constant a d as: 



(29) 



a d = Pd.S/2 



but every positive smaller value is fine; we remind that a d is approximatively the absolute value of some 
interactions on <9i?22, so in general a d is a small positive number. 
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Given a set of edges E let us define He as the contribution to the Hamiltonian on the edges E; 
precisely we have: 



(30) 



For the next lemma sec also (H],|9). 



H E {cr) = - J x , y o- 



{x,y}£E 



Lemma 4.4. Let 5 G (0, 1/2). There exists /3i(<5) such that for all [3 > (3\{8) the following upper bound 
is satisfied 



(31) 



7rp({\M(dB 2U a)\ < l-S})+7r ({\M(dB 2{ i +1) ,a)\ < 1 - 6}) < e^ 1 "' 1 



where K x = d2 d ~ l 5 2 p 2 d . 



Proof. If the absolute value of the magnetization is less than (1 — 5) on the set dB 2 i or on the 
set <9-B 2 (i+i) then, using Lemma 4.2 and Lemma 4.3, there exists a Peierls surface 7 where at least a 



percentage Spd of edges are unsatisfied; moreover all these edges have interactions that are close to 1. 
The satisfied edges are all on 7 n (dbB 2 i U dbB^) and the absolute value of the interactions of these edges 
is small (the values are close to or sad if e € dbB 2 i and close to zero if e £ dbBa). For a fixed Peierls 
surface 7 if a configuration a respects the previous conditions about the satisfied and unsatisfied edges 
then we write a ~ 7. To bound the probability we consider a spin flip on the interior of 7. This changes 
the Hamiltonian only on the edges of 7 (see fig. |]) leaving the same values on all the other edges; this 
gives us an easy bound on the probability of all the configurations a ~ 7 (see 0, §). On the edges of 

+ , + 



+ 



+ 



- + 
+ + 



+ 



+ 



+ 



+ 



+ - + 



+ 



+ 



+ 



Figure 2. The spin-flip of the region inside the circuit 7 



7 = 7(0") -where 7 is the Peierls surface of Lemma 4.2 and 4.3- the value of the Hamiltonian has these 
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bounds: 

(32) 5p d \<y\ - a d \j\ -l<H 7 (a)< Sp\j\ + a d | 7 | + 1, 

for all cr verifing \M (dB2i, cr)| < 1 — 6 or M(<9£? 2 (;+i), cr) | < 1 — 5. The inequalities are clear because 
all the edges in 7 \ (dbB2i U dbB^i) -they are at least <5p|7|- are unsatisfied; in the first inequality we 
are considering that all the other edges in 7 n (9fci?2/ U dbB^) are satisfied so they can give a contribute 
to the Hamiltonian that can be only larger than — ad| 7 | — 1; vice versa in the second inequality we are 
considering that also the edges in 7 n (dbB 2 i U dbB^) are unsatisfied. Analogously, taking the spin-flip of 
the configuration on the interior of 7, we obtain 

(33) -6p d \<y\ - a«l|7l - 1 < H^a) < -5p\~/\ + a d \j\ + 1, 
So 

(34) H(a) - H(f y a) = H^(a) - # 7 (/ 7 <r) > 2Sp d \ 1 \ - 2a d | 7 | - 2. 
Now, using relation (p9|), one has 

(35) H(a)-H(f^a)>dp d \j\-2. 

Let us fix a Peierls surface 7 such that \j \ (c^i^i U dbBn)\ > <5prf It | ■ For all the a having all the edges 
on 7 \ (dbB2i U dbBn) unsatisfied we find this bound (we are using inequality (|35|)) 

(36) J ^ = ^^<^ (5pM - 2) - 

7173 (/ 7 er) e -P H (M 

For a fixed 7 we perform the sum over all the a ~ 7 (i.e. all the configurations a having r )\{dbB2i UdbB^i) 
as a subset of boundary of clusters): 



(37) M°) < e-W^-V £ ^(f 7 a) < e 



/3(5p d | T |-2) 



since the function / 7 : Q — ► f2 is injective, and so the sum Ylcr~-y 7T 0if / y (7 ) ^ s a probability (therefore a 
number belonging to [0, 1]). Now, taking all the possible 7 which fulfill the conditions of Lemma L2J-4.3, 
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and using inequality ([37]), we obtain: 

7T P ({\M(dB 2l ,a)\ < 1 -8}) +wp({\M(dB 2{l+1) ,a)\ < 1 - 6}) 



(38) 



7:\7\>Sp d \dbB 2 i\ n>5p d \d b B 2 i\ 



where Kq is an uniform upper bound for the number of the Peierls surface of size n in any dimension d 
(see 0, [221). Now, if (3 is large enough, we can bound the ( |38| ) with the expression 



exp (~/35 2 p 2 d \dB 2l \/2) = exp (~/3S' pp^'d^ 1 + p(l)), 



2 „2 d-l j/id-l 



where p(l) is a o(l d 1 ). So for (3 large enough we can choose K\ = d2 d 1 S 2 p 2 l . rj 



Remark 4.5. // the Markov chain starts with the stationary measure then also the event (S' ; + (4(5) U 



S t (45)) c has probability smaller than e Kl @ l at every time ( e N (we are also using Lemma J^A_ to 
connect the magnetization with the events S^(A5) and 5 ; _ (4<5)J. So we have: 



(39) 



where K Y =K 1 /A = d2 d - A 6 2 pj 



^((5+wu5 r wr)< e -^ d 



Let us define the following function $ : Q — » D, inside a cube Bu with interactions verifying the 
conditions 1, 2 and 3 of the second section: 

-a(Rv) if v g i?2i 

(40) = <( CT (i? w ) if v g (B 4 , \ £ 2( ) 

ct(u) if u ^ B 4 ; 

and R is the reflection with respect to the center of the cube B41. It is easy to see that $ 2 = 1 (the 
identity function). By construction we have that, for every S g (0,1), if a g S + (5) then 6 5' _ (^) 

and if <r g S~(6) then $(cr) g S + (5). The function <!> has these important properties: 



(41) 



M(dB 2l ,${a)) = -M(dB 2h a) and M(dB 2{l+1) , $(a)) = M{dB 2(l+1)) a) 
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Also every spin inside the cube B41 sees the same next neighbored except the spins on the boundaries 
dB2i, <9£?2(i+i) and dB^. We notice that an edge e £ dbB 2 i is satisfied for the configuration a if and only 
if Re is satisfied for the configuration $(cr); this observation will be essential in the proof of the next 
lemma. 

Lemma 4.6. Let's 8 e (0, 1/2). Then there exists s E (0, 1 — 25) and fctf) such that for all j3 > /3 2 (5) .- 
(42) a) lim^(S+(<S)) = l, 

I — >oo 

^ 6) hmsup 7^1^r = 0, 

;-oo Trp{S l (5))e K2 P l 

where the positive constant K 2 is equal to d2 d {\ — s — 28)8pd- 



Proof. By Lemma 4.4 and Remark 4.5 we know that the probability of S, .(5) U S l (8) goes to one 



when I goes to infinity. Now we will prove that the measure will be concentrated on the set S^(8) for / 
tending to infinity. 

Using the function $ : Q — > il we want to compare the probability of any configuration a £ Sf(S) 
with the probability of the configuration $(cr) 6 Sf(8). For every realization of the interactions that 
respects our constraints (see Section 2 points 1-3) we have an uncertainty about the Gibbs factor of order 
e ± ' 3C ', where C is a constant independent of I, because the contribution to the Hamiltonian inside the 
cube B41 depends on the fact that the interactions are values belonging to a small interval of length l/l d 
(the inverse of the volume so they are inversely proportional to the number of interactions inside the cube 
B41). We do not control the value of the interactions on the boundary of the faces (they have dimension 
(d — 2) so the number of interactions are 0(l d ~ 2 )). Therefore, including the two kinds of mistakes, we 
can calculate the Hamiltonian with an uncertainty of order l d ~ 2 . Let us consider the configuration $(cr) 
where a is any configuration in S^(8). As noted, just before this lemma, we have that every satisfied 
(resp. unsatisfied) edge e S dbBii for the configuration a corresponds to a satisfied (resp. unsatisfied) 
edge Re G dbB 2 i for the configuration $(c). Moreover, if e £ dbB 2 i has J e > (resp. J e < 0), then 
Re has Jn e < (resp. J Re > 0), this is a consequence of the definition of the interactions inside the 
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cube B<u. Now we give some rigorous upper and lower bounds for the expression 7T | a($(cr))/7r l a(cr) when 
a G Sf(5); we will use some other bounds for the the Hamiltonian on the edge set d b B 2 i. Doing algebraic 
calculations we find: 



(44) \\-\-*) °*\hBn\ < H 9bB2l (a) < (±- I 



+ ss) a d \d b B 2l \ VtreSf(5,t); 



where the lower (upper) bound is obtained with a configuration a having K(d b B 2 i, cr) = 1 — 6 with all 
the positive (negative) edges e = {i,j} = d b B 2 i having OiGj = 1. We also find: 

(45) (~ + I - s5j a d \d b B 2l \ < H dbB2l {<j) < (-1 + | + sj a d \d b B 2l \ Vrr 6 S+{5). 

So, if the parameter 6 is small, all the allowed values of Hg b B 2l on Sf (5) (resp. on S^(S)) belong to 
a small interval. Using the bounds (f5|)-(ff5|) we can find these useful bounds for the difference of the 
Hamiltonian: 

(46) (-l + s-2sS)a d \d b B 2 i\ < H^b^i) - H db B 2l ((X 2 ) < (-1 + s + 25)a d \d b B 2t \ 

for every o \ S Sf(6), a 2 G S^(5). Now we choose s € (0,1) and 6 G (0,1/2) such that (— 1 + s — 2s6) and 
(— 1 + s + 26) will be negative. In the end, see Theorem |4.8| , we will prove that the different conditions on 
the parameters 6 and s can be simultaneously satisfied. One has \d b B 2 i\ = C^^ 1 — d2 d l d ~ 1 ; so, using 
inequality ( jiff ) we find 

(47) e P(-l+s-2s8)a d C d l^ < < e P(-l+s+28)a d C d l^ for ^ ff £ g+ {§ y 



So we deduce from the first inequality of (47) that 



(48) hmsup^» = 0, 



this means that for large j3 -where we can apply Lemma 4.4 and Remark fLq- we obtain 



(49) lim tt p (S+(6)) = 1. 
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Using the second inequality in ( f47| ) we deduce that 

where K 2 can be chosen as any positive constant larger than d2 d (l — s — 28)ad — d2 d ~ 1 (l — s — 25)pdS, 
we take K 2 — d2 d (l — s — 2S)pdS. This concludes the proof, rj 

Lemma 4.7. Let's 5 G (0, 1/2). There exists ^(5) such that for all (3 > /3s(5) the probability that either 
(51) <r{t) G S+{6) and a{t + 1) e Sf(<J) 



(52) cr(t) e Sf(5) and a(t + 1) G S+{5) 
is smaller than e" K3 ^ 1 where K3 — spdSd2 d ~ 5 . 

Proof. The first step SW1) in the evolution of configuration of edge r\ is really simple because all the 
variables {r] e } ee E A are independent. The probability to have on an edge e = {v, w} € E\ rj e (t + 1) = 1 
depends only on the spins a v (t) and a w (t). More precisely if J v ,w°~v(t)o~w(t) > then rj v ^ w takes value 1 
with probability 1 — e -P\ J v,vj\- jf J v w a v (t)a w (t) < r] VtW is equal to with probability 1. 

Let us consider the case {a(t) G Sf(5)} — > {a(t + 1) G 5 ; + (6)}. Since cr(t) is in Sf(S) there are at 
least (1 — 5/2)|9bi?2i| edges e = on d b B 2 i having a v a w = —1. Among these edges at least 

(l-S/2)\d b B 2l \ - \d b B 2l \/2 

are negative edges. So at least (1 — 5)\df,B2i\/2 are negative edges having probability 1 — p e to be occupied 
(i.e. rj e = 1). There is the following uniform lower bound for the probability that one of these edges is 
occupied: 

(53) p e > inf{l - e _/3|Je| : e G <9 6j B 2 z is a negative edge } > 1 - e^ 80 ". 
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Now we use the standard large deviation estimation for a Bernoulli distribution; let us suppose that 
are i.i.d. random variable with P p [X\ = 0) = 1 — p and P p {X\ = 1) = p then 



(54) P p X i < xN "j < e- I(x ' p > N if x < E(Xi) = p, 
where the rate function I(x,p) = xlog(^) + (1 — x) log(^5^). For any fixed x G (0, 1) 

(55) Urn -^4^ r = 1. 

V ' P-i- (l-x)log(l-p) 

So for any fixed x G (0, 1 ) there is p G (0, f ) large enough such that 

(56) I(x,p) > i-^log ' 



2 1 -p 

Now we can estimate the probability that in the step SWf) of SW dynamics there are less than 



|i^|9bi?2/| occupied edges when at least ^-\dbB 2 i\ edges can be occupied 



\e£A 



1 1 — 5 ■ » ' ' ' -f3sa d /2^- 5 



(VT) P„ \ > , v. < g — I^BmIJ < exp(-J(-,l - e -f«-^)__|^B 2 ,|) J 

we are using ([si]) with p that is any positive number less or equal to inf{l— e _/3 Jc ' : e G <9(,£?2i is a negative edg 
we have chosen p = 1 — e~P sa <*/ 2 ; A C d\,B 2 i is the set of edges that can be occupied and \A\ > ^-^-\dbB>2i\. 
Now, for a fixed S G (0, 1/2) and for (3 large enough, using (Bq) we can bound the probability in (|57j) with 



(58) P p I > > e < L-i^Bail ) < cxp(-sa d (3\d b B 2l \/32). 



We notice that if there are at least ^^|9f,i?2i| occupied negative edges in the step SW1) then in the 
second step SW2) for every occupied negative edge e = {v,w} the product a v a w will be equal to —1, 
so K{dbB2i, cr(i)), see definition (pd|), will be less than (1 + 5)/2 and the configuration a(t + 1) is not in 
S^(S) when 5 G (0, 1/2). In this case the constant K$ can be chosen smaller or equal to sadd2 d ~ 4 (see 
inequality (|58|)). 

The other case {<j(t) G S^^)} — > {<r{t + 1) G ^(J)} is analogous and in both transitions we can take 



if 3 = sadd2 d ~ 4 — spd5d2 d ~ 5 (we have used (|29|) for the last equality), rj 

In the following theorem we will denote with L the length of the side of the box A = and l(L) 
will denote an opportune function such that one can find, for L large enough, inside Aj, a box B 4 i^ 
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verifying all the hypotheses on the interactions written in the second section. It is clear that a strictly 
increasing function l(L) exists because of finding the right interactions is an event of probability larger 
than zero and the interactions are i.i.d. random variables. In the next theorem we denote with the 
box Z d n [— L,L] d and (3(6) — max{/3i(<$), 02(6), /?3(<5)}. In the next theorem we consider j3 larger than 
/3(d) so are verified simultaneously all the inequalities of previous lemmas. 

Theorem 4.8. Let 6 e (0, 1/2), (3 > 0(6) , and f L ,s : Q,a l — » R 



(59) 



f L ,s(u)) = l s - ( a)H 



then, the spectral radius Rl goes to one when L increase to infinity Qj-almost surely and so 



(60) 



lim T exp . L — oo Qj - a.s. 

L — >oo 



Proof. First we need to show that it is possible to select the constant K\, K 2 and K3 of Lemmas 4.4- 
4.7 and Remark |4.5| such that they satisfy the inequalities < K2 < K\ and K% < K3. We show that 
there are some values of the parameters 5 and s such that all the inequalities between K\, K% and K3 
are verified. The inequalities to be verified are: 

1 - s - 28 > (=>- Kx > 0) 
S 2 p 2 d > 16(1 - s - 26) (=* Ki > K 2 ) 
( 61 ) { sp d 8 > 64(1 - s - 25) K 3 > K 2 ) 

< s < 1 
< <5 < 1/2 



We remind that is a fixed geometrical constant in the open interval (0, 1) (see Lemma 4.2 and 
Lemma |4.3| ); it is immediate to see that for every 5 G (0,1/2) there are some values of the parame- 
ter s such that all the inequalities are verified (a non empty interval of values) . One can choose any value 
of s in the open interval 
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s < 1 — 2(5 from the first inequality in (|6l|) , s > max | Pd |^ 64 (1 ~ 2<5), f — 2(5 — 4^ | from the second and 
third inequalities in (|6l|). It is immediate to see that also the other inequalities are verified and 

max( r 64 — (1 - 25), 1 - 26 - ^\ < 1 - 26 
\p d <5 + 64 v ; ' 16 J 

for every <5 £ (0, 1/2). In the rest of the proof we will omit the dependence of the eigenvalues from the 
parameter 5. 

Now we use the Dirichlet's formula (see equation (||)) and the variational characterization of Ai to 
prove that the spectral gap goes to zero when L increase to infinity. The usual variational characterization 
of the second eigenvalue gives: 

(62) Ax = 1 - inf j f J, y : < Var w (f) < oo 



Var„(f) 



and also 



(63) A* = 1 - inf \ -^41 : < VoivCf) < oo 



Var„(f) 

for every <eN; where t is the number in which the configuration is updated with SW algorithm. 
If we can prove that there exists an increasing sequence of times = t(l(L)) j oo such that 

(64) ton °' Lj f f L ] -0 Qj-a.8. 

then we get that the spectral gap goes to zero Qj-a.s. (see for similar arguments Q). 

Let us define /; = Jl = ^-g~(L)> then for every time t G N we have the following upper bound 



(65) Cf /(t) < tTra(n \ {Sr U 5,+)) + t max M| w (cj, J) + * max M§ w {u, J) 



where we are supposing that there is a box Bn with all the interactions described in the second section. 
The first summand in RHS of (^5|) is an upper bound for the probability that in the time {1, . . . , t} the 
chain visits a configuration out of U Sf ; this is a consequence of the initial measure that is equal to 
the stationary measure 7rg, so at any time the associated measure of the Markov chain is equal to the 
upper bound is a consequence of sub-additivity. The second summand in RHS of ([35]) is an upper bound 
of a transition from a state of to a state of . Analogously the third summand in (p5|) control the 
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probability of a transition from to . We are imposing that every time that there is a transition 
in s G {l,...,t} from S[~ to Sj or a transition from 5 ; + to ST" or the chain visits a configuration in 
51 \ U S*^) then the function (l(Sf)(uit) — l(Sf)(uJo)) 2 of formula (||) is equal to 1 where lo s is the 
configuration of the process at time s G {l,...,t} (this is an obvious upper bound). To find an upper 



bound for the RHS of (|65|) we will use Lemma 4.4 and Lemma 4.7, so one obtains that 



(66) C fl j l (t)<t[e- K ^ ld - 1 + 2e- 



For C/, ,/,(()) = Var(f), using Lemma 4.6, we obtain the following lower bound 



C M M>\e K ^ 



So 



(67) < 2t[e-^- K ^ t + 2e-^- K ^ ld ~\ 

c fi,fA Q ) 

where the constant (K\ — K2) and (K3 — K2) are positive. So one can choose an increasing sequence of 
times ti, depending on I, such that ti -^>i^oo 00 and 

Jim = 0. 

C/ (>/! (0) 

To finish the proof we stress that one can find a sequence of boxes (with side length L) such that 
just a finite number of these boxes have not inside a square B41 with all the interactions verifying 
hypotheses 1-3 of the second section. So Rl goes to one Qj-a.s. and the exponential autocorrelation 
time goes to infinity using the relation (pi]), rj 

Some remarks: We can consider a limit dynamics when the temperature parameter T goes to zero; 
it is clear that in this case there is a different behavior for the ferromagnetic systems with respect to 
the frustrated ones. For the ferromagnetic systems all the configurations that are not ground states for 
the Hamiltonian are transient states and the unique invariant measure of the associated Markov chain is 
uniform on both the ground states (all the spins equal to 1 or -1). Moreover the ground states form an 
irreducible class so with probability one the Markov chain visits all these ground states (the states are 
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finite because we are considering just a finite box A). For frustrated systems there are different invariant 
measures and the ground states are divided in several irreducible classes; so it is immediate to see that in 
this case the spectral gap is equal to zero because there are states that are not reachable, starting from 
some states, in spite of a positive probability for some invariant measures. 

We want to stress that not only our proof does not work for the ferromagnetic systems but also 
we believe that really there is a different behavior for all the small values of the temperature between 
ferromagnetic and frustrated systems. Our conjecture is that for every ferromagnetic system there exists 
To such that for each T £ [0, To) the spectral gap of the associated Markov chain is uniformly bounded 
away from zero independently of the size of the systems (the size of the box A) . This is not in contrast 
with because a Markov process on an infinite space of states can be non-ergodic also if all the terms 
in the sequence of its approximations on finite spaces of states have uniform bounds (away from zero) for 
the spectral gap. 

Vice versa in this paper we have just proved that for every frustrated disordered system, where the 
interactions have £ supp Qj, the spectral gap shrinks to zero when the size of the system increases to 
infinity and the temperature is small enough. So we believe that there is a different behavior for small 
values of the temperature between the ferromagnetic and the disordered not ferromagnetic systems. 

Following one could prove that for small and large values of f3 there is a different behavior of the 
convergence to the equilibrium with respect to the distance D p . The distance D p between two measures 
/ii and /12 defined on { — 1, 1} A is: 

(68) D p (fn, fj, 2 ) = sup ^2 ^2 IM^A+i) - H2(<?A+i)\ P W 

l£A A,i£AcA<r A e{-l,l} A 

where A is a connected vertex set; p £ (0, 1) is a parameter. For small values of (3 and of the parameter 
p we obtain that 



(69) 



D p ( Mo Mf w *, M A,j,/3) <Me~ 
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where c > can be chosen independent of A. For large values of (3 the relation (^9|) is not true for any 
fixed value of p. 

5. Conclusions 

We have proved some bounds on the spectral radius of the Swendsen-Wang dynamics and we found 
a different behavior for large or small values of the temperature parameter. Our upper bound on the 
spectral gap is not explicit but if one needs an explicit evaluation one should only perform our calculation 



in Theorem 4.8 preserving all the coefficients. Also one could find, for every distribution of the interactions 



J , the optimal values of the parameters 8 and s to estimate, using the Dirichlct's formula, the spectral 
gap of Ml J. 

A comparison with the SW dynamics of the ferromagnetic Ising model would be interesting. We believe 
that the ferromagnetic systems carry out better performances than the disordered not ferromagnetic one 
at every value of the temperature 1//?; moreover we think that the spectral gap in the ferromagnetic case 
does not shrink to zero when the size increases to infinity if the temperature is small enough. 

Acknowledgments. / wish to thank Charles Newman and Alan Sokal for beneficial conversations and 
for encouraging me to work on this paper. I thank Jeff Cheeger which have given me some ideas to prove 
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